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Heavy Quark Potentials and Quarkonia Binding 
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Abstract. I review recent progress in studying in-medium modification of inter-quark forces at finite tem- 
perature in lattice QCD. Some applications to the problem of quarkonium binding in potential models is 
also discussed. 
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1 Introduction 

The study of in-medium modifications of inter-quark forces 
at high temperatures is important for detailed theoretical 
understanding of the properties of Quark Gluon Plasma 
as well to detect its formation in heavy ion collisions. In 
particular, it was suggested by Matsui and Satz that color 
screening at high temperature will result in dissolution 
of quarkonium state and the corresponding quarkonium 
suppression could be a signal of Quark Gluon Plasma for- 
mation 1 . 

Usually the problem of in-medium modification of inter- 
quark forces is studied in terms of so-called finite tem- 
perature heavy quark potentials, which are, in fact, the 
differences in the free energies of the system with static 
quark anti-quark pair and the same system without static 
charges. Alternatively this problem can be studied in terms 
of finite temperature quarkonium spectral functions [2 
which were also discussed during this conference by 
Karsch, Hatsuda and Petrov |5j . Recently substantial prog- 
ress has been made in studying the free energy of static 
quark anti-quark pair which I am going to review in the 
present paper. An important question is what can we learn 
about the quarkonium properties from the free energy of 
static charges which will be discussed at the end of the 
paper. 



2 The free energy of static charges 

Following McLerran and Svetitsky the partition function 
of the system with static quark anti-quark (QQ) pair at 
finite temperature T can be written as 



ZQQir,T) = {W{r)WHO))Z{T), 



(1) 



with Z(T) being the partition function of the system with- 
out static charges and 



W{x) = V exp{ig / drAo(r, a; 

^0 



is the temporal Wilson line. L{x) — TrM^(a;) is also known 
as Polyakov loop and in the case of pure gauge theory it 
is an order parameter of the dcconfinement transition. As 
the QQ pair can be either in color singlet or octet state 
one should separate these irreducible contributions to the 
partition function. This can be done using the projection 
operators Pi and Pg onto color singlet and octet states 
introduced in Refs. [311]. Applying Pi and Ps to Zggir, T) 
we get the following expression for the singlet and octet 
free energies of the static QQ pair 



exp(-Fi(r,T)/r) 



1 TvPiZnnir^T) 



Z{T) TrPi 
\T:T{W{r)W\Q)) 



(3) 



exp(-F8(r,r)/r) = 



1 TvP^ZQQjr.T) 
Z{T) 



TtPr 



)) 



= \(TrW{r)TvW\Q)) - ^Tr{W{r)W\Q)). (4) 
8 24 

Although usually Fi g, is referred to as the free energy 
of the static QQ pair, it is important to keep in mind 
that it refers to the difference between the free energy of 
the system with static quark anti-quark pair and the free 
energy of the system without static charges. 

As W{x) is a not gauge invariant operator we have to 
fix a gauge in order to define Fi and Fg, . As we want that 
Fi and Fg have a meaningful zero temperature limit we 
better to fix the Coulomb gauge because in this gauge a 
transfer matrix can be defined and the free energy differ- 
ence can be related to the interaction energy of a static 
QQ pair at zero temperature (T = 0). Another possibil- 
ity discussed in Ref. [5] is to replace the Wilson line by 
a gauge invariant Wilson line using the eigenvector of the 
spatial covariant Laplacian ^ . For the singlet free energy 
both methods were tested and they were shown to give 
numerically indistinguishable results, which in the zero 
temperature limit arc the same as the canonical results ob- 
tained from Wilson loops. The interpretation of the color 
octet free energy at small temperatures is less obvious and 
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Fig. 1. The color singlet free energy in quenched |1()I28| (top) 
and three flavor |14| (bottom) QCD. The solid black line is the 
T — singlet potential. 



will be discussed separately. One can also define the color 
averaged free energy 



eM-Favir, T)/T) 



1 Tr(Pi+P8)^go(r,r) 
Z{T) Tr(Pi + Fg) 



= i(TrTy(r)TrVFt(0)), 



(5) 



which, is expressed entirely in terms of gauge invariant 
Polyakov loops. This is the reason why it was extensively 
studied on lattice during the last two decades. The color 
averaged free energy is a thermal average over the free 
energies in color singlet and color octet states 



exp(-F,„(r,T)/T) - 
i exp(-Fi(r, T)/T) + | exp(-F8(r, T)/T). 



(6) 



Therefore it gives less direct information about medium 
modification of inter-quark forces. Given the partition func- 
tion Zqq (r, T) we can calculate not only the free energy 
but also the entropy as well as the internal energy of the 
static charges 



dF,{r, T) 



dT 



'^^^'^^ dT \ Z{T) 

= F,{r, T)+TS,{r, T), 
i = 1, 8, av. 



, (7) 



(8) 
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Fig. 2. The effective screening radius versus T /Tc I14j . 

3 Free energy of a static QQ pair and 
screening of inter-quark forces at high 
temperatures 

Perturbatively the quark anti-quark potential can be re- 
lated to the scattering amplitude corresponding to one 
giuon exchange and in the non-relativistic limit it is given 

by 



V{r) = {T''T'')g^ 



(9) 



Here Doo{k) is the temporal part of the Coulomb gauge 
gluon propagator and in general it has the form 

Dooik) = (k2 + 77oo(k))-i. 

Furthermore the averaging over color gives {T°-T^) — —4/3 
for the color singlet and {T°-T^) = +1/6 for the color octet 
case. At zero temperature the polarization operator TTqo 
gives rise only to running of the coupling constant g — g{r) 
(recall that as = (7^/(47r)). But at finite temperature T 
it has a non-trivial infrared limit noQ{k — > 0) = rrv^ = 
gTy/Nc/3 + Nf/6. Therefore at distances r > 1/T the 
potential has the form 



V{r, T) = (r°T'')f — exp(-m£,r 



47rr 



(10) 



The singlet and octet free energies defined in the pre- 
vious section can also be easily calculated in leading order 



perturbation theory. Again because of i7oo(fc = 0) = mjj 
one has 

Pi3(r,T) = (-^ i)^exp(-m^r). (11) 



At leading order the singlet free energy has exactly the 
same form as the potential and has no entropy contribu- 
tion. This is the reason why the free energies of static QQ 
pair were (mis)interpreted as potentials. At next to lead- 
ing order which is 0{g'^) the free energies have the form 



4 1. 



9 



exp(— mjjr) 



(12) 



3 ' 6 47rr ^ ^ ^ ' 37r 

and the entropy contribution —TS appears (recall Eq.[7I) 
For the singlet case the entropy has the form 



Siir,T) 



3ttT 



(1 — exp(— TO£ir)). 



(13) 
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Fig. 3. The color octet free energy in quenched QCD |23| . The 
solid black line is the zero temperature singlet potential. 
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It has the asymptotic value of ^-^f- at large distances and 
vanishes for rT <C 1. Similarly one can calculate the color 
octet entropy to be 

^8(r,r)-^(l + iexp(-mz,r)). (14) 

Contrary to the color singlet case it does not vanish at 
small distances. We can also calculate the internal energy 
which for color singlet state, for example, can be written 

as 

Ui{r,T} = --- — exp(-m£,r) exp(-m£,r), 

d 47rr Sir 

(15) 

and unlike the free energy vanishes at large distances (at 
least to order g^). Using Eqs. Hll|) one can easily get 
the perturbative result for the color averaged free energy. 
For rT > 1 the exponentials in Eq. ^ can be expanded 
and we arrive at the well known leading order result jHllHI 

1 g4 

Fav{r, T) = ^ exp(-2mz3r). (16) 



4 Numerical results on the free energy of 
static QQ 

4.1 Color singlet free energy 

In this section I am going to review recent lattice results on 
the free energy of a static quark anti-quark pair. I will start 
the discussion with the case of the color singlet channel. 
The color singlet free energy has been extensively stud- 
ied only during the last three years. Presently results are 
available for SU(2) and SU(3) gauge theories pminnrra 
[T5| as well as in two flavor ^2] and three flavor QCD 
[Tl] . While for pure gauge theories these studies are very 
systematic and lattice artifacts are under control, for full 
QCD they are still in the exploratory stage. 

In Fig. n the singlet free energy for SU(3) gauge the- 
ory (QCD without dynamical quarks or quenched QCD) 
is shown for different temperatures together with the zero 
temperature quark anti-quark potential. For temperatures 
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Fig. 4. The ratio of color singlet and color octet free energies 

below the transition temperature Tc ~ 270MeV the free 
energy rises linearly with the distance r signaling confine- 
ment. Above deconfinement T > Tc the free energy has 
a finite value at infinite separation indicating screening. 
One can also see from the figure that at short distances, 
rT <^ 1, the free energy is temperature independent in 
the entire temperature range and equal to the zero tem- 
perature potential. The singlet free energy for three flavor 
QCD is also shown in Fig. ^ The main difference com- 
pared to the case of SU(3) gauge theory is the fact that 
the free energy reaches a constant value at all tempera- 
tures. At low temperatures this is interpreted as string 
breaking, the flux tube breaks if enough energy is accu- 
mulated to create a light quark anti-quark pair which with 
the static QQ could form a static-light meson, i.e. when 
V{r = rscr) = E^heavy-ught ( ^66 e.g. discussiou in Ref. [El 
) . The distance rscr , where the free energy effectively flat- 
tens off depends on the temperature, it becomes smaller 
as the temperature increases. Therefore it is interpreted 
as an effective screening radius and is shown in Fig. |21 
At low temperatures, T < Tc it has a value of about 0.9 
fm and rapidly decreases near the transition point. While 
close to Tc the effect of dynamical quarks is important for 
the value of the screening radius, at high temperatures 
the value of the screening radius is similar in quenched 
and full QCD. 

4.2 Color octet free energy 

The color octet free energy is shown in Fig.Ofor quenched 
QCD. At short distances it is repulsive as expected from 
perturbation theory. Above the deconfinement tempera- 
ture it has strong temperature dependence which presum- 
ably comes from the entropy contribution and is present 
even at short distances. At high temperatures this is also 
expected from perturbation theory (see previous section). 
Above deconfinement the color octet free energy has the 
same large distance asymptotic value as the color singlet 
free energy Fs{r ^ oo,T) = Fi{r ^ oo,T) = F^{T). 
This is intuitively expected, at large distances the quark 
and anti-quark are screened by their respective "clouds" 
and do not know anything about their relative color ori- 
entation. One should note that also below Tc the octet 
and singlet free energies become equal at large distances ( 
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Fig. 5. Tlie iiybrid potentials in quenched QCD from Ref. 1191 . 
^9^('") — ^i''') is the singlet potential and the dashed lines are 
predictions from string model. 



compare Figs. ^ and O though there is no particular phys- 
ical reason for this. I will discuss this problem at the end 
of this subsection. These features of the color octet free 
energy are present also in the case of three flavor QCD 

m 

Perturbation theory predicts that at high tempera- 
tures we expect for the ratio of 6Fi_g,{r, T) = {Fi^si^, ^) ~ 
i^oo(T)) the following dFi{r,T)/ Fs{r,T) - -8. In Fig. H 
the lattice data in SU(3) gauge theory are confronted with 
these expectations. As one can see from the figure the data 
for this ratio are close to —8 for temperatures T > 2Tc. 

While at high temperatures the meaning of the no- 
tion of color octet free energy is clear, the meaning of this 
quantity at low temperatures ("confinement region") is 
less evident. To understand the problem let us first dis- 
cuss the spectrum of static quark anti-quark free energies 
at zero temperature. The lowest energy level of a static QQ 
pair is the one where the quark and anti-quark are color 
singlet state. The corresponding energy as a function of 
the quark anti-quark separation is the singlet static poten- 
tial or simply the static potential determined in terms of 
Wilson loop and used extensively in potential models (see 
Refs. |17lll8n. There are also higher energy levels, whose 
energy functions are called hybrid potentials for which the 
gluon fields between the static charges are in excited state 
(or in other words the string formed between the quark 
and anti-quark is excited) .VT^ld ■,20.- The spectrum of hy- 
brid potentials is shown in Fig. [3 The hybrid potentials 
are labeled by the angular momentum projection of the 
gluon field configurations on the quark anti-quark axis, 
L = 0,1,2 (denoted as U, 77, A), CP (even, g, or odd, 
u) and the reflections properties with respect to the plane 
passing through the quark anti-quark axis (even, -I- , or odd, 
-) • The most distinct feature of the hybrid potentials 
is the different slope at small distances, i.e. in contrast to 
the singlet potential the hybrid potentials are repulsive. 
They have a shape which is similar to the shape of the 
octet free energy at low temperatures shown in Fig. O It 
has been shown that at short distances hybrid potentials 
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Fig. 6. The quenched running coupling constant at finite tem- 
perature in the color octet case (23; . The thick black line rep- 
resents the lattice data at zero temperature on the running 
coupling. The thin black line is the running coupling derived 
from Coulomb plus linear parametrization of the zero temper- 
ature potential. Finally the thin gray line is the 3-loop running 
coupling in qq scheme |24| . 

correspond to the perturbative color octet potential |21j . 
This means that at short distances hybrid potentials cor- 
respond to a state in which the static QQ pair is in octet 
state and the net color charge is compensated by soft glu- 
ons field which make the whole object color singlet (ob- 
viously only singlet objects can exist in the confinement 
region). At very small distances the soft gluon field is de- 
coupled and the energy is dominantly determined by the 
large repulsive interaction of the static quark and anti- 
quark 

The correlators which enters the definition of the color 
singlet and octet free energy have the following spectral 
representation 

(TrM^(r)Trl^t(o)} = J^,, e-^-f"--^)/^ (17) 

{TrWir)W^O)) = E„ c„e-^"('-'^)/^, (18) 

where En denotes the different energy levels of the quark 
anti-quark system: singlet potential, hybrid potentials, sin- 
glet potentials plus glueballs, etc. The weights c„ in gen- 
eral are different from one and may have non-trivial r- 
dependence Because of asymptotic freedom ci should 
approach unity at short distances, while c„>i should van- 
ish; at short distances perturbation theory can be applied 
and the correlator {TiW {r)W'' (0)) gives the singlet po- 
tential. The tendency of ci approaching unity is clearly 
seen in the lattice data presented in Ref. [22]. Since the 
gap between the zero temperature (ground state) poten- 
tial and the lowest hybrid potential is large for not too 
large distances the color octet free energy is expected to 
be given by 

,-Fsir,T)/T^l_^^-mr)/T^ E,(r)^V{r), (19) 

i.e. it is determined by the singlet potential V{r). This is 
the reason why at very large distances the color singlet and 
color octet free energies are equal in the low temperature 
region (compare Fig.^and Fig.O. On the other hand for 
sufficiently small distances 1 — ci (r) ~ and the octet free 
energy is expected to be determined by the lowest lying 
hybrid potential Fs{r,T) ~ £'2(r) -f rin8. 
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Fig. 7. Tlie internal energy (to) and entropy (bottom) of static 
QQ pair at infinite separation in quenched |28| . 2 flavor |15| 
and 3 flavor QCD 14 . The filled and open black circles is the 
free energy in 2 flavor and quenched QCD respectively. 




0.2 0.4 0.6 0, 



Fig. 8. The internal energy in quenched QCD il2,. The solid 
black line is the T = potential. 



out that aeffi — a^ffs- In Fig. Elthe running coupling 
constant is shown for the octet case (the resuhs in the 
singlet case are essentially identical). At short distances 
the effective coupling constant is temperature indepen- 
dent and coincides with the zero temperature result. At 
larger distances its deviates from the zero temperature re- 
sult and after approaching a maximum it drops because of 
the onset of screening. Note that for temperatures close to 
but already above Tc the effective running coupling con- 
stant follows the zero temperature running coupling up 
to distances of about 0.3 fm. At such distances the run- 
ning coupling is not controlled by the perturbation theory 
but rather by the linear part of the potential which gives 
the nearly quadratic rise of this quantity. Thus some non- 
perturbative confining physics survives deconfinement. 



5 Entropy and internal energy of a static 
quark anti-quark pair 

As it was noticed in section 2, given the partition func- 
tion Zqq we can calculate the entropy and internal energy 
difference of the system with static charges and the same 
system without them, which for the sake of simplicity is 
called entropy and internal energy of QQ. Using Eqs. Q 
and © the entropy and internal energy has been calcu- 
lated for infinite separation in quenched two flavor 
PI and three flavor UH QCD. The results are shown in 
Fig. [7| Both the entropy and the internal energy show a 
very large increase near Tc- For quenched QCD the inter- 
nal energy has been calculated for any separation r [T^ 
and the results are shown in Fig. |S1 One can see that with 
the exception of the small distance region where the inter- 
nal energy coincides with the zero temperature potential 
the internal energy is larger than the free energy. In the 
high temperature limit the internal energy has no constant 
piece at large distances proportional to the temperature, 
therefore it is tempting to interpret the free energies as 
potentials. However, the large increase of U{r,T) near Tc 
make such interpretation problematic. 



4.3 Running coupling constant at finite temperature 

To study how the free energies approach the zero temper- 
ature limit as well as to make contact with perturbation 
theory at short distances it is convenient to introduce the 
effective running coupling constant ae//(r). This quan- 
tity can be also used to quantify the strength of interac- 
tion at least within the perturbative framework. At zero 
temperature the most convenient way to introduce the ef- 
fective running coupling constant is through the force be- 
tween quark and anti-quark, ae//(r) = (3/4) • r^{dV/dr). 
This definition avoids many problems in the perturba- 
tive calculation of the effective coupling constant f2[;. 
One can define similar quantities at finite temperatures 
Pll^: aeff i,8(r,T) = (3/4,-6) •r2(dFi,8(r,T)/dr) For 
the color singlet case in quenched QCD ae// was dis- 
cussed in detail in Ref. where it was also pointed 



6 Quarkonium binding at finite temperature 

Following the suggestion by Matsui and Satz P the prob- 
lem of quarkonium binding at finite temperature has been 
studied using potential models with some phenomenolog- 
ical screened potential (see e.g. 26,27,28 ) which led to 
the conclusions that J/ip dissolves in the Quark Gluon 
Plasma at temperatures close to Tc- More recently the free 
energy has been used as the potential in the Schroedinger 
equation and it was found that J/ip can survive only to 
temperatures l.lTc As the free energy contains an 
r-dependent entropy contribution the validity of this ap- 
proach is doubtful. Finally, very recently the internal en- 
ergy calculated in Ref. ^21 was used as a potential in 
Schroedinger equation [301 and it was found that J/ tp can 
survive till l.TTc which is not inconsistent with lattice cal- 
culations of J/ip spectral function PIE]- However, to test 
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the validity of potential models it is not sufficient to make 
statements about the dissolution temperature of a given 
quarkonium state, one should investigate the change in 
the properties of heavy quarkonium bound states. For ex- 
ample, potential models with screening predict a decrease 
in the quarkonium masses. The most convenient way to 
compare the prediction of potential models with direct 
calculation of quarkonium spectral functions is to calcu- 
late the Euclidean meson correlator at finite temperature 
(see contribution by Mocsy to this proceedings ^TB). This 
quantity can be reliably calculated on the lattice. The ba- 
sic idea is to use the model spectral function 

i 

containing bound states (resonances) and continuum |31[ 
1^2 . For a given screened potential one can solve the 
Schroedinger equation and determine the radial wave func- 
tion (or its derivative) at the origin i?i(0) and the binding 
energy Ei. The parameters of the model spectral func- 
tions can be related to these quantities, FiiT) ^ |_R(0)|^, 
Mi = 2mc,b -I- Ei and acquire temperature dependence 
because of the temperature dependence of the potential. 
Here rric^b is the constituent quark mass of c- and b-quarks. 
The threshold of the continuum so{T) can be related to 
the asymptotic value of the potential at infinite distance 
So(T') = 2mc,b -I- Vao{T). Having specified completely the 
spectral function of the model the Euclidean correlator 
can be calculated 



/■CSC 

G(t,T) = / duja{uj,T) 



cosh(w(T - i/(2r))) 

sinh(cj/(2T)) 



(20) 



and compared with lattice results [35] ■ The Euclidean cor- 
relator can be reliably calculated on lattice while extract- 
ing the spectral function from it is quite difficult. Using 
this approach some of the features of the quarkonium 
correlators, e.g. the enhancement of the scalar correlator 
above deconfinement temperature observed in lattice cal- 
culation, can be understood. More detailed discussion on 
this topic is given in Ref. |32j . 



7 Conclusions 

Free energies of static QQ pairs have been extensively 
studied on lattice and provide useful tool to study in- 
medium modification of inter-quark forces. Many body ef- 
fects exert large influence on the free energy thus making 
the simple picture, where the temperature dependence of 
the free energy reflects the screening of the two-body po- 
tential, not applicable. Many body effects are most promi- 
nent close to the transition temperature. For temperatures 
not too close to the transition temperature free energies of 
static quark anti-quark pairs could provide useful qualita- 
tive (though not quantitative) insights into the problem of 
quarkonium binding in Quark Gluon Plasma. It is inter- 
esting to note in this respect that a new model approach 
to the problem of heavy quark potentials at flnite temper- 
ature was recently proposed in Ref. . 



Acknowledgements 

This work was partly supported by U.S. Department of 
Energy under contract DE-AC02-98CH10886. P.P is a Gold- 
haber and RIKEN-BNL Fellow. The authour would like 
to thank F. Zantow for correspondence and A. Patkos 
for careful reading of the manuscript and valuable sug- 
gestions. 

References 

1. T. Matsui, H. Satz, Phys. Lett. B 178, 416 (1986). 

2. T. Umeda, K. Nomura, H. Matsufuru, 
|arXiv:hep-lat /0211003| 

3. M. Asakawa, '1'. Hatsuda, Phys. Rev. Lett. 92, 012001 
(2004) arXiv:hep-lat/0308034 . 

4. S. Datta, F. Karsch, P. Petreczky, I. Wetzorke, Phys. Rev. 
D 69, 094507 (2004) 

5. Contributions by F. Karsch, T. Hatsuda and K. Petrov to 
these proceedings. 

6. L.D. McLerran, B. Svetitsky, Phys. Rev. D 24, 450, (1981) 

7. L.S. Brown, W.I. Weisberger, Phys. Rev. D 20, 3239, 
(1979) 

8. S. Nadkarni, Phys. Rev. D 33, 3738 (1986) 

9. O. Philipsen, Phys. Lett. B 535, 138 (2002) 

10. O. Kaczmarek, F. Karsch, P. Petreczky, F. Zantow, Phys. 
Lett. B 543, 41 (2002) 

11. S. Digal, S. Fortunato, P. Petreczky, Phys. Rev. D 68, 
034008 (2003) 

12. O. Kaczmarek, F. Karsch, P. Petreczky, F. Zantow, Nucl. 
Phys. Proc. Suppl. 129, 560 (2004) 

13. O. Kaczmarek, F. Karsch, F. Zantow, P. Petreczky, Phys. 
Rev. D 70, 074505 (2004) 

14. P. Petreczky, K. Petrov, Phys. Rev. D 70, 054503 (2004) 

15. O. Kaczmarek, F. Zantow, progress; also O. Kaczmarek, 
these proceedings 

16. S. Digal, P. Petreczky, H. Satz, Phys. Lett. B 514, 57 
(2001) 

17. G. S. BaU, Phys. Rept. 343, 1 (2001) 

18. N. Brambilla et al, hep-ph/0412158 

19. C. J. Morningstar, K. J. Juge, J. Kuti, Nucl. Phys. Proc. 
Suppl. 73, 590 (1999) 

20. K. J. Juge, J. Kuti, C. Morningstar, Phys. Rev. Lett. 90, 
161601 (2003) 

21. N. Brambilla, A. Pineda, J. Soto, A. Vairo, Nucl. Phys. B 
566, 275 (2000) arXiv:hep-ph/9907240 . 

22. O. Jahn, O. Philipsen, Phys. Rev. D 70, 074504 (2004) 

23. O. Kaczmarek, F. Karsch, P. Petreczky, F. Zantow, work 
in progress 

24. S. Necco, R. Sommer, Phys. Lett. B 523, 135 (2001); Y. 
Sumino, Phys. Rev. D 65 054003 (2002) 

25. F. Zantow, these proceedings 

26. F. Karsch, M. T. Mehr, H. Satz, Z. Phys. C 37, 617 (1988). 

27. G. Ropke, D. Blaschke, H. Schulz, Phys. Rev. D 38, 3589 
(1988) 

28. T. Hashimoto et al., Z. Phys. C 38, 251 (1988) 

29. S. Digal, P. Petreczky, H. Satz, Phys. Rev. D 64, 094015 
(2001) 

30. E. V. Shuryak, I. Zahed, Phys. Rev. D 70, 054507 (2004); 
C. Y. L. Wong, hep-ph/0 408020. 

31. E. V. Shuryak, Rev. Mod. Phys. 65, 1 (1 993). 

32. A. Mocsy, P. Petreczky, hep-ph/04 lT262l 

33. Yu. A. Simonov, hep-ph/0502078. 



